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Abstract. The curl-modes of Cosmic Microwave Background (CMB) 
polarization probe horizon-scale primordial gravitational waves related 
to inflation. A significant source of confusion is expected from a leas¬ 
ing conversion of polarization related to density perturbations to the 
curl mode, during the propagation of photons through the large scale 
structure. Either high resolution CMB anisotropy observations or 21 cm 
fluctuations at redshifts 30 and higher can be used to delens polariza¬ 
tion data and to separate gravitational-wave polarization signature from 
that of cosmic-shear related signal. Separations based on proposed lens- 
ing reconstruction techniques for reasonable future experiments allow the 
possibility to probe inflationary energy scales down to 10^^ GeV. Beyond 
CMB polarization, at frequencies between 0.01 Hz to 1 Hz, space-based 
laser interferometers can also be used to probe the inflationary gravita¬ 
tional wave background. The confusion here is related to the removal of 
merging neutron star binaries at cosmological distances. Given the low 
merger rate and the rapid evolution of the gravitational wave frequency 
across this band, reliable removal techniques can be constructed. We 
discuss issues related to joint constraints that can be placed on the infla¬ 
tionary models based on CMB polarization information and space-based 
interferometers such as the Big Bang Observer. 


1. CMB: At Present 

The cosmic microwave background (CMB) is now a well known probe of the early 
universe. The temperature fluctuations in the CMB, especially the so-called 
acoustic peaks in the angular power spectrum of CMB anisotropies, capture the 
physics of primordial photon-baryon fluid undergoing oscillations in the potential 
wells of dark matter (Hu et al. 1997). The associated physics — involving the 
evolution of a single photon-baryon fluid under Compton scattering and gravity 
— are both simple and linear, and many aspects of it have been discussed in 
the literature since the early 1970s (Peebles &: Yu 1970; Sunyaev &: Zel’dovich 
1970). The gravitational redshift contribution at large angular scales (e.g., Sachs 
&: Wolfe 1968) and the photon-diffusion damping at small angular scales (e.g.. 
Silk 1968) complete this description. 

By now, the structure of the first few acoustic peaks is well studied with 
NASA’s Wilkinson Microwave Anisotropy Probe (WMAP) mission (e.g., Spergel 
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et al. 2003), while in the long term, ESA’s Planck surveyor^, will extend this 
to a multipole of ~ 2000. Given the improved resolution and better frequency 
coverage, a variety of studies related to secondary anisotropies, e.g., Sunyaev- 
Zel’dovich effect (Sunyaev & Zel’dovich 1980; Cooray, Hu & Tegmark 2000) 
is expected. Beyond temperature fluctuations, detections of the polarization 
anisotropy, related to density or scalar perturbations, have now been made by 
DASI (Kovac et al. 2003; Leitch et al. 2004) and CBI (Readhead et al. 2004; 
Cartwright et al. 2005), while the temperature-polarization power spectrum is 
measured by WMAP (Kogut et al. 2003). 

An interesting result from WMAP data is related to the very large ( > 10°) 
angular scale polarization signal, which probes the local universe and associated 
astrophysics instead of physics at the last scattering surface probed with polar¬ 
ization at tens of arcminute scales. This local universe contribution arises when 
the universe reionizes again and the temperature quadrupole begins to rescatter 
to produce a new contribution to the polarization anisotropy at angular scales 
corresponding to the horizon at the new scattering surface (Zaldarriaga 1997). 
The large-scale excess signal measured in the WMAP temperature-polarization 
cross-correlation power spectrum is interpreted as rescattering with an optical 
depth to electron scattering of 0.17 ± 0.04 (Kogut et al. 2003). 

Such an optical depth suggests early reionization; for example, if the uni¬ 
verse reionized completely and instantaneously at some redshift, the measured 
optical depth suggests a reionization redshift of ~ 17 (±5). Such a high optical 
depth, in the presence of a > 1% neutral Hydrogen fraction from z ~ 6 Sloan 
quasars (Fan et al. 2002), suggests a complex reionization history. A complex, 
and patchy, reionization, however, is expected if the reionization is dominated 
by the UV light from first luminous objects, though to explain the high level of 
the WMAP’s optical depth requires high star formation efficiency at redshifts of 
order 15 with a possibility for a population of metal-free massive Pop HI stars 
(e.g., Cen 2003). 

In general, one can only extract limited information related to the reion¬ 
ization history from polarization peak at tens of degrees or more angular scales 
(e.g., Hu & Holder 2003). One reason for this is the large cosmic variance associ¬ 
ated with measurements at multipoles of ~ 10. To extract detailed information 
of the reionization process, one can move to temperature fluctuations generated 
during the reionization era at angular scales corresponding to few arcminutes 
and below. At these small angular scales, potentially interesting contributions 
from high redshifts include a scattering contribution related to moving electrons 
in ionized patches and fluctuations in the moving electron population (e.g., San¬ 
tos et al. 2003), and a thermal Sunyaev-Zel’dovich effect related to the first 
supernovae bubbles (Oh et al. 2003). The same redshift ranges are expected to 
be illuminated in the near-infrared band between 1 to 3 ^m in the form of spatial 
fluctuations to the cosmic IR background (e.g., Santos et al. 2002; Salvaterra 
et al. 2003; Cooray et al. 2004). Potentially interesting studies to extract more 
information related to the reionization process include brightness temperature 
fluctuations of the redshifted 21 cm line emission by neutral HI prior to and 
during Hydrogen (e.g., Zaldarriaga et al. 2004), and cross-correlation studies 


^ http://astro.estec.esa.nl/Planck/ 
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Figure 1. Power spectrum for the temperature anisotropies in the 
fiducial ACDM model with r = 0.1. In the case of temperature, the 
curves show the local universe contributions to CMB due to gravity 
(ISW and lensing) and scattering (Doppler, SZ effects, patchy reion¬ 
ization). We refer the reader to Cooray, Baumann & Sigurdson (2004) 
for a recent review on large scale structure contributions to tempera¬ 
ture anisotropies. 


between small angular scale effects in CMB and maps at other wavelengths 
(Cooray 2004a; Cooray 2004b; Salvaterra et al. 2005). Upcoming experiments 
such as the South Pole Telescope and the Atacama Cosmology Telescope, com¬ 
bined with near-IR and 21 cm maps, will provide some of the first possibilities 
in this direction. 


2. Inflation: CMB and Space-Based Laser Interferometers 

We hrst summarize basic aspects related to inflation under the assumption of a 
single scalar field model. Our derivations, however, are applicable to multi-field 
theories with few alterations (Lidsey et al. 1997). In general, constraints on infla¬ 
tion are obtained by comparing theoretically calculated scalar and tensor metric 
— gravitational waves — perturbations to direct (in the case of space-based 
detectors) and indirect (such as CMB) measurements of these perturbations. 

The tensor and scalar perturbations are determined by the evolution of 
the scalar held, (/>, (the ’inflaton’ field) and the scale factor, a, with time. The 
equation of motion for a scalar field with a self-interaction potential, V{4>), is 
given by, 

'(j) + + U' = 0, 


( 1 ) 
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Figure 2. CMB polarization anisotropies in the gradient (E) mode 
due to scalars and curl (B) modes due to the inflationary gravitational 
wave background (curve labeled ’IGB’, with the number in parenthesis 
showing the normalization in terms of the tensor-to-scalar ratio, which 
in terms of the energy scale of inflation is 2 x 10^® GeV (0.1), top curve, 
and 2x 10^® GeV (10“^), bottom curve). Note the effect due to reioniza¬ 
tion where rescattering produces new large angular scale anisotropies. 
For optical depths at the level of 0.1, the peak of the power related to 
IGW curl modes is at tens of degree angular scales and not at the de¬ 
gree scale {I ~ 100) corresponding to the bump at recombination. Note 
that this enhancement, amounting to over two magnitudes in power at 
I ~ few is significant and aid potential detectability of the IGWs via 
polarization. The IGW detection is confused with cosmic shear conver¬ 
sion of a fractional E-mode —> B-mode by the intervening large scale 
structure. The residual noise curve related to cosmic shear is after do¬ 
ing lens-cleaning in polarization maps with an arcminute scale CMB 
experiment with instrumental polarization sensitivity at the level of 
1 ;uK -y/sec. For reference, we show an estimated level of contamina¬ 
tion from dust and synchrotron emission (at 150 GHz); this noise level 
can be improved with foreground cleaning in multifrequency data with 
sensitivity better than the foreground noise contributions. The curve 
labeled Planck is the noise curve related to upcoming Planck (HFI) 
polarization observations as a function of each multipole. 


where H is the Hubble parameter and the prime refers to differentiation with 
respect to the field, (j)- The evolution of the Hubble parameter is determined 
by the energy density, p, which we suppose is dominated by the energy density 
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associated with the inflaton field, 




Stt Stt 



( 2 ) 


where Mpi = = 1.2 x 10^® GeV. The initial curvature redshifts away by the 

time the scales corresponding to the current Hubble horizon left the inflationary 
horizon. In general, an exact analytic solution for the spectrum of scalar and 
tensor perturbations does not exist. In order to facilitate analytic solutions, one 
generally resorts to the ’slow-roll’ expansion, where the equations of motion are 
expanded in terms of the inflaton potential and its derivatives (Liddle, Parsons 
& Barrow 1994). Motivation for this expansion can be understood by noting 
that the equation of motion for the inflaton field is formally equivalent to the 
equation of motion of a mass moving under a potential V with damping due to 
the cosmological expansion. The slow roll expansion supposes that the motion 
of the field reaches a ’terminal velocity’ in which the potential energy dominates 
over the kinetic energy and the acceleration of the field is negligible. In order to 
quantify these assumptions, the slow-roll parameters are used 


^ IbTT V H j ’ 


and 




^ 4 ^ 

Stt V ' 


(3) 


Under the assumption that e, \ri\ <C 1, the approximate ’slow-roll’ evolution 
equations are 




Stt 


r, 


with 3H ^ + V' 


0 . 


(4) 


On a given length scale, scalar perturbations start as quantum fluctuations in 
otherwise classical fields well within the horizon. Here, flat space-time quantum 
field theory is appropriate. As the scale factor quasi-exponentially increases, 
the length scale is stretched beyond the inflationary horizon. The quantum 
fluctuations are then amplified and become classical perturbations to the metric 
(Guth & Pi 1982). In order to study scalar and tensor perturbations, we consider 
their power spectra. It can be shown that within the slow roll approximation 
the amplitude of power spectra can be written as (Stewart & Lyth 1993, Lidsey 
et al. 1997), 



1 

1287r 


Ps{k) = 

TT m 2 e ^ 

SM^i U'2 ’ 

(5) 


16 9 

128 V 


Prik) = 

O H 

7rM2i 

3 VpV 

(6) 


where S and T refer to scalar and tensor perturbations respectively. The right 
hand sides of the above equations are to be evaluated when a given comoving 
length scale crosses the inflationary horizon with /cq = aH. In order to quan¬ 
tify how these power spectra change with scale, one considers their power-law 
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spectral indices (Liddle and Lyth 1992), 


ns{k) 

nrik) 


dlnPs 

dink 


1 — 6e + 2 ri, 


din Pt 
dink 


(7) 

( 8 ) 


If the potential is flat enough, higher order spectral parameters, which are higher 
order in the slow roll parametersm, can be ignored (Stewart & Lyth 1993). As 
long as the slow roll parameters can be considered constant, we are able to 
write the power spectra as power-laws with spectral indices given by the above 
expressions. It can then be shown that, 

Ae^ Ar]^ 0{e^)Alnk. (9) 


Currently, CMB data probe up to Imax ~ 1000 corresponding to Aln/c ~ I. 
Current observations constrain e < 0.1 at scales corresponding to the current 
Hubble horizon. Therefore, although CMB observations do cover a few orders of 
magnitude in wavenumber, to within second order in the slow roll parameters, 
they can be considered constant over the range probed by the CMB. This implies 
that the perturbation power spectra will appear as approximate power laws, 
expanded around some ’pivot’ wavenumber, ko (Leach et al. 2002), 

Ps{k) ^ Psiko) , and P^ik) ^ Priko) . (10) 

In the case of WMAP these parameters are evaluated at ko = 0.002Mpc“^ 
(Spergel et al. 2003; Peiris et al. 2003). Using the expressions for the amplitudes 
and spectral indices in terms of the slow roll parameters, we are able to solve 
for the potential and its first two derivatives and thereby use observations to 
constrain the form that this potential can take. We have, 


V S 

- Mi. 

128 P 

(11) 

|U'| p 


(12) 

V" p 

- ■^7V(ns' - 1 - 3nT)Mpi, 

(13) 

Ps 

^ —WriT, 

(14) 


and all of the spectral parameters are to be evaluated at the pivot wavenumber, 
ko- The last expression is known as the consistency relation and is true only in 
single-held inhationary model. The expression relaxes to an inequality in the 
case of multi-held models. It is a consequence of the fact that both scalar and 
tensor spectra are determined from the dynamics of a single scalar held, so they 
are not independent (Lidsey et al. 1997, Liddle & Lyth 2000). 
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2.1. CMB 

Though acoustic oscillations in the temperature anisotropies of CMB suggest 
an inflationary origin for primordial perturbations, such as in the WMAP data 
(Peiris et al. 2003), it has been argued for a while that the smoking-gun signature 
for inflation would be the detection of a stochastic background of gravitational 
waves (e.g., Kamionkowski & Kosowsky 1999). Since CMB probes horizon to 
super-horizon scales, inflation is the only known mechanism to casually pro¬ 
duce a stochastic background of gravitational waves at such large wavelengths 
(e.g., Starobinsky 1979). Note that the amplitude of the gravitational wave 
background is highly unknown. Since the amplitude of these inflationary grav¬ 
itational waves is proportional to V, the value of the inflaton potential V{4>) 
during inflation, the amplitude of gravitational-wave induced B-mode polariza¬ 
tion anisotropies directly constrains the energy scale of inflation Relative 

to COBE normalization, the relation between the energy scale and the ratio of 
tensor-to-scalar fluctuations is = 3.0 x 10“^ (T/S')^/^Mp; (Turner & White 
1996), The current limit from CMB anisotropy data, based on large angular- 
scale CMB temperature fluctuations, suggests an upper limit of order 2 x 10^® 
GeV (Melchiorri & Odman 2003). 

Instead of temperature fluctuations (where dominant scalar perturbations 
confuse the detection of tensor contribution), the gravitational-waves can be 
studied through their distinct signature in the CMB polarization in the form of 
a contribution to the curl, or magnetic-like, component (Kamionkowski et al. 
1997; Seljak & Zaldarriaga 1997). There is no contribution from the dominant 
scalar, or density-perturbation, contribution to these curl modes in linear theory. 
Secondary effects related to density perturbations, however, do produce signals 
in the curl-mode and confuse the detection of inflationary gravitational-wave 
signal. 

Among these, the few percent conversion of E to B from cosmic shear mod¬ 
ification to the polarization pattern by the intervening large scale structure is 
important (Zaldarriaga & Seljak 1998). These lensing-induced curl modes intro¬ 
duce a noise from which gravitational waves must be distinguished in the CMB 
polarization. 

Gravitational leasing of the CMB photons by the mass fluctuations in the 
large-scale structure is now well understood (e.g., Hu 2000; Hu & Cooray 2001). 
The leasing effect can be described through a remapping process involving an¬ 
gular deflections resulting along the photon path; 

0(h) = 0(h-|-V(()) 

= 0(h) + Vi((>(h)V*0(h) + iVi</>(h)Vj(/>(h)V*V%(h) + ... , (15) 

where the deflection angle on the sky is given by the angular gradient of the 
leasing potential, a(h) = V(/>(h), which itself is a projection of the gravitational 
potential, (see e.g. Kaiser 1992), 

(/>(m) = —2 / dr <h(r, ihr). (16) 

Jo dA{r)dA{ro) 

The quantities here are the conformal distance or lookback time, from the ob¬ 
server (r) and the analogous angular diameter distance dA- The leasing potential 
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in equation (16) can be related to the well known convergence generally encoun¬ 
tered in conventional leasing studies involving galaxy shear (Kaiser 1992) 


K(m) 


^V2(/>(rh) 



^ jA{r)dA{ro 

dAiro) 


-V^d>(r, rhr), 


(17) 


where note that the 2D Laplacian operating on <1> is a spatial and not an angular 
Laplacian. Expanding the leasing potential to Fourier moments, 

r ^2i 

<t>A) = J . (18) 

we can write the usually familiar quantities of convergence and shear components 
of weak leasing as 


-n(n)±trAti) = (19) 

Taking the Fourier transform, as appropriate for a flat-sky, we write 

0 (li) = J dh 0 (h)e-*'i-" 

= 0(ii)-y ^0(ii')^(ii,iio, (20) 

where 

1 r d^" 

L(ii,iiO = 0(1, -1/) (h -1,0.1/+- y ^^(1/) (21) 

x 0 *(l," + l,'-l,)(l,".l,O(li" + li'-li)-li'. 


The leased power spectrum, according to the present formulation, is dis¬ 
cussed in Hu (2000) and we can write 


Cf = 


1 - 


dH 




C? + 


dH 


1 


- iii ■ ■. 


ii'i 


( 22 ) 


Since the leasing effect involves the angular gradient of CMB photons and 
leaves the surface brightness unaffected, its signatures are at the second order in 
temperature. Effectively, leasing smoothes the acoustic peak structure at large 
angular scales and moves photons to small scales. Note that the approach above 
treats the deflection angle as a perturbative parameter. Such an approximation 
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breaks down at angular scales corresponding to roughly the rms deflection angle 
of a CMB photon (~ a few arcminutes). While the approach can be extended 
to higher order (e.g., Cooray 2004c), the numerical calculation becomes increas¬ 
ingly time consuming. The development of numerically fast, but still accurate, 
approaches (e.g., Challinor & Lewis 2005) will aid in analysis data from future 
high resolution CMB experiments. To put rough estimates on the expected 
changes due to lensing, the mean gradient is of order ~ 15 ;uK arcmin“^ and 
with deflection angles or order 0.5 arcmin or so from massive clusters, the lensing 
effect results in temperature fluctuations of order ~ 5 to 10 ^K. In Fig. 1, we 
show the effect of lensing on temperature fluctuations where power is transferred 
from acoustic peaks to the damping tail of anisotropy spectrum. 

To calculate the lensing contribution to the polarization, we follow the no¬ 
tation in Hu (2000). Introducing, 

±X(fi) = ±X(fi + V</.) (23) 

±X(h) + Vi0(h)V'±X(h) + iVi<()(h)V,</.(fi)V*V^±X(h), 

where = Q dz iU represent combinations of the Stokes parameters. The 
lensing effect is again to move the photon propagation directions on the sky. In 
Fourier space, assuming flat-sky, we can consider the E- and B-mode decom¬ 
position introduced in Kamionkowski et al. (1997; Seljak & Zaldarriaga 1997) 
such that ±-^(1) = E(l) ± iB{l) and 

±X(1) = ±X(l)-y ^±X(li)e±2*(‘^-i-‘^-)L(l,li). (24) 

Since primordial H(l), due to the gravitational wave contribution, is small, we 
can make the useful approximation that 

E{1) = E(l)-y'^[E(li)cos2(^^-(^i)]</.(l-li)(l-li).li 

m = -y ^[E(li)sin2((^i,-<^i)]</>(l-li)(l-li)-li, (25) 

Under this approximation, the lensed polarization power spectra can be ex¬ 
pressed in terms of Cf^ and, in the case of B-modes, we can write 




r<BB,IGW 



(27r)2^|i-iil 


[(l-ll)-ll]2Qf^[l 


cos(4(^iJ] , 


(26) 


BB IGW 

where (7^ ’ is the contribution due to IGWs and the contribution in the 

second line is the confusion that must be removed. To achieve this separation, 
one must obtain information on the integrated potential out to last scattering 
surface. One can effectively extract information related to the lensing potential, 
and the integrated dark matter density field responsible for the lensing effect, 
via quadratic statistics applied to CMB temperature and polarization maps (Hu 
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Figure 3. CMB B-mode polarization. The curve labeled ‘IGWs’ is 
the IGW contribution with a tensor-to-scalar ratio of 0.1 with (solid 
line; r = 0.17) and without (dashed line) reionization. The curve la¬ 
beled ‘lensing’ is the total lensing confusion to B-modes. Thin lines 
show the residual B-mode lensing contamination for removal of with 
lensing out to Zg while thick lines show estimates of the residual confu¬ 
sion from GMB experiments alone using quadratic estimators with an 
ideal noise-free experiment (dashed line), and likelihood methods using 
a high resolution/sensitivity experiment (dot-dashed). The noise curve 
of latter this experiment, with 2 arcminute beams and a pixel noise 
of 0.25 /rK-arcminute, is the curve labeled ‘High-Res’. Bias-limited 
lensing information to Zg > 10 improves upon the limit to the IGW 
amplitude based on quadratic statistics, and the likelihood level can 
be reached with Zg ~ 30. If Zg ~ 100 an additional order-of-magnitude 
in the IGW amplitude could be probed. The curve labeled ‘Low-Res’ 
is the noise curve for a lower resolution CMB polarization experiment 
with 30 arcminute beams sufficiently sensitive to detect IGWs when 
paired with a cosmic 21-cm lensing reconstruction. See, Sigurdson &: 
Cooray (2005) for details. 


2001a, Hu 2001b; Hu & Okamoto 2003; Cooray &: Kesden 2003; Kesden et al 
2003) or likelihood-based techniques (Hirata &: Seljak 2003). 

In figure 2, we show the contribution from dominant scalar modes to the 
polarization in the gradient component and from gravitational-waves to the curl 
polarization. Unlike the B modes generated by tensor perturbations (the IGWs), 
is dominated by larger amplitude scalar perturbations. The expected 
few-percent conversion of E modes creates a large signal in the B-mode power 
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spectrum (Zaldarriaga & Seljak 1998). For tensor-to-scalar ratios T/S below 
2.6 X 10“^ or below 4.6 x 10^® GeV the IGW signal is completely confused 
by the lensing contaminant (Lewis et al. 2002). To bypass this limit one must 
separate the lensing induced B-modes from those due to IGWs (Kesden et al. 
2002; Knox & Song 2002; Seljak & Hirata 2004). Glearly, the lensing confusion 
could be exactly removed if one knew the three-dimensional distribution of mass 
out to the GMB last-scattering surface. However, knowledge of the projected 
quantity 4){n) is sufficient. One way to estimate <?!)(h) is by using quadratic 
estimators or maximum likelihood methods to statistically infer the deflection- 
angle field given temperature and polarization anisotropies (Hu 2001b; Hu & 
Okamoto 2002). These method allow one to detection down 10^® GeV (Kesden 
et al. 2002; Knox & Song 2002). In the case of likelihood methods, one does not 
have an obvious ultimate limit, however, the residual lensing noise is still above 
the detector noise even for highest sensitive polarization maps. While this is 
encouraging, it is useful to consider alternative techniques to reduce the confu¬ 
sion. A useful way to estimate ())(h), but out to a high-redshift, is by observing 
the weak-lensing distortions of the shapes of objects of a known average shape 
at redshift Zg- As the GMB source redshift is at 1100, observations of the weak 
lensing of galaxies, out to a redshift of ~ 1-2, can not be used to effectively 
delens GMB maps because a large fraction of the lensing contamination (55% 
at I = 1000) comes from structure at z > 3. A possibility exists in the form 
anisotropies in the 21 cm brightness temperature fluctuations during the era and 
prior to reionization (Sigurdson & Gooray 2005). 

For example, using a quadratic reconstruction of the deflection field of 21 cm 
anisotropies out to a source redshift of 30, in conjunction with a low-resolution 
but high signal-to-noise GMB polarization experiment, could detect T/S > 2.5 x 
10“^ or > 2.6 X 10^® GeV if 21 cm fluctuations measurements are limited 
to I ~ 5000. If the resolution in 21 cm observations is improved so that one 
can measure 21 cm brightness temperature anisotropy spectrum out to I 10^ 
one can improve down to T/S > 1.0 x 10“® or > 1.1 x 10^® GeV. This 
limit is comparable to using GMB lensing information alone from high sensitive 
maps at 2 arcminute resolution (Seljak & Hirata 2004). To achieve such a 
resolution one requires roughly a space mission with an aperture of diameter 
~ 6(100 GHz/i/) meters. In the case of GMB lens-cleaning with information 
from 21 cm observations, as the IGW signal is at multipoles of 100 (in the 
case of recombination bump) and/or 20 (in the case of reionization with optical 
depths around 0.1), the GMB data need not be high resolution. Since lensing 
information is not extracted from GMB observations, the optimal observing 
strategy is to integrate over a few square degree patch of the sky as proposed in 
Kamionkowski & Kosowsky (1998). Thus, it is conceivable that the Inflationary 
Probe of NASA’s Beyond Einstein Program can be designed in combination with 
a cosmic 21-cm radiation experiment. Such a cosmic 21-cm radiation experiment 
could be extremely exciting for reasons involving 21 cm fluctuations alone where 
one probes the small scale structure of the Universe at high redshifts. For a bias- 
limited reconstruction out to a ~ 100, the limit on the tensor-to-scalar ratio 
is 7.0 X 10“® or > 6.0 x 10^^ GeV. For Zg ~ 200, the maximum redshift 
where 21-cm fluctuations are expected to be nonzero, one could probe down to 

> 3 X 10^^ GeV. 
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2.2. Space-based Laser Interferometers 

Beyond CMB polarization data, there is long-term interest, as part of NASA’s 
Beyond Einstein Program, to launch a space-based laser interferometer (Big 
Bang Observer) to directly detect the presence of relic stochastic gravitational 
background today from inflation. The CMB detection and the detection of the 
laser interferometers provide unique opportunities to study properties related to 
inflationary physics. To understand requirements for the direction detection, we 
summarize some of the key aspects related to the stochastic background. 

The evolution of tensor perturbations, for scales inside the horizon subse¬ 
quent to inflation, is governed by the massless Klein-Gordon equation 


d^hk 


+ 2 


1 da dhk 
a dr dr 


+ k^hk — 0 


(27) 


subject to the boundary conditions, hk{0) = Prik)^^"^ and hk{0) = 0. r is 
conformal time, which is related to cosmic time by the differential relation, 
dt/dr = a. We define gk{T) = hk{T)a{T). With this definition, we are able to 
rewrite the above equation, 


d^Qk 

dr"^ 



1 d^a\ 


0 


(28) 


There are two limiting behaviors for gk‘. before horizon entry gk oc a ^ hk = 
constant; after horizon entry gk will oscillate which implies hk will oscillate with 
an amplitude decreasing as 1/a. Therefore, the current spectrum of gravitational 
waves is determined by the primordial power spectrum and by the rate at which 
scales enter the horizon (i.e. the evolution of the scale factor); h = hoak/a. At 
horizon crossing, note that k = akHk- 

During radiation domination, H oc a“^, so that k = l/ak and during matter 
domination H oc so that k = . From these relations, we find that, 

hk oc k~^ (Radiation Domination) (29) 

hk oc k~‘^ (Matter Domination) (30) 

In order to quantify these scalings, a transfer function is used (Turner, White 
& Lidsey 1993). We note that during the matter dominated epoch, a oc and 
the above equation becomes a Bessel equation (see, discussions in Pritchard & 
Kamionkowski 2004). Requiring that the solution approach unity at r = 0, we 
find that 

hkir) = 3:^i^ (31) 

Using this, we define the transfer function, T{k/keq) (where keq is the wavenum¬ 
ber of matter-radiation equality ~ 6.22/i^ x 10“^ Mpc“^) as. 


{hkikrY) 


((3ii(fer)/(A:r))2) 


T{k)^ 


(32) 
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where the average is taken over several periods around the time that we want to 
evaluate the transfer function. From our scaling arguments, the transfer function 
must take the form of 


T{k) = ^Co + Cikjkeq + C2ik/keqy 


(33) 


The spectrum as seen by a gravitational wave detector is (Turner, White & 
Lidsey 1993), 


flcwik) 



Tjkjke,? 

{kroY 


(34) 


where tq ^ k, 10^ Mpc for Hq = 70 km/s Mpc. flow is the ratio of 

the energy density in gravitational energy to that of the closure density, pc ~ 
10“2® gm/cm^. In order to determine the values for ci,C 2 , and C 3 we need to 
numerically integrate the evolution equation and determine T[k). Since space- 
based gravitational wave detectors that are under consideration for detection 
of inflationary gravitational waves have sensitivities around a frequency of ~ 
O.OlHz (~ 10^^ Mpc“^), our analysis considers only wave numbers in excess of 
keq- Therefore, the values determined by Turner, White and Lidsey (1993; cq = 
l,ci = 1.34, C 2 = 2.5) are applicable, even though they assumed the universe 
was filled with only matter and radiation. 

Here, we have ignored one important detail; the right hand side of time 
evolution equation is not strictly zero. Instead, tensor perturbations are sourced 
by anisotropic stress in the cosmic fluid. The largest source of such stresses comes 
from the free-streaming of relic neutrinos (Weinberg 2004). It can be shown 
that these neutrinos cause the amplitude of gravitational waves that enter the 
horizon well before matter radiation equality to be damped by a factor of ~ 0 . 8 . 
This will influence the amplitude of gravitational waves detectable with space- 
based detectors. Similarly, the amplitude seen by CMB polarization is damped 
between a few to a few ten percent with longest waves damped less. 

In Fig. 4, we show the amplitude of tensor fluctuations related to a back¬ 
ground of inflationary gravitational waves as a function of the gravitational wave 
frequency. We plot roughly 25 orders of magnitude in the frequency from super¬ 
horizon to horizon scales probed by CMB to a million-km wavelength modes 
detectable with space-based interferometers. For reference, we also show vari¬ 
ous existing limits on the stochastic backgrounds based on calculations related 
to the Big Bang Nucleosynthesis (Allen 1996) to timing-limit for millisecond 
pulsars (e.g., Thorsett & Dewey 1996). An extended discussion related to such 
limits are in Seto & Cooray (2005). The left-panel of Fig. 3 shows Hgw which is 
the fractional density contribution from a background of stochastic gravitational 
waves with the density pcw- 




1 dpGW 
Pcdlogf ■ 


(35) 


For reference, in the right-panel, we also plot the dimensional strain hdf) given 
by the h‘l{f) = fSh{f) where Sh{f) is the power spectral density of gravita¬ 
tional waves, {h{f)h{f')) = Sh{f). The two quantities, Dgw and are related 
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through (Maggiore 2000), 


^Gwif) = . (36) 

In addition to single detector noise levels, in Fig. 4, we also show the im¬ 
provement when the signal from two independent detectors can be combined 
through a correlation analysis between the two detectors (e.g., Maggiore 2000). 
Such a correlation increases sensitivity to a stochastic background such that, 

= 1.12 X 10-S(/)MI (^)''‘ (37) 

where A/ is the bandwidth over which the signals can be correlated, T is the 
integration time, and F is a filling factor, which accounts for the fact that a pri¬ 
mordial stochastic background will be isotropic on the sky, but the detector will 
only be sensitive to a fraction of the sky. For omni-directional interferometers, 
such as the Laser Interferometer Space Antenna (LISA), F = ‘Ijh, while hn{k) 
is the char acteristic strain due to detector noise. The single-detector sensitivity 
is hn{k)/F. Here, y{f) is the overlap-reduction factor that takes in to account 
the relative orientation and configuration of two interferometer systems. This 
correction results in an overall reduction of noise from the naive estimates that 
do not take in to account detail configuration and motion of the detectors with 
respect to each other in the orbit (e.g., Ungarelli & Vecchio 2001). 

For a correlation analysis, the increase in sensitivity is under the assumption 
that the detector noises are independent between the two detectors, while the 
only correlation expected is due to stochastic signals such as inflation. We 
note here that for the single detector case, the minimum observable strain is 
independent of the integration time while for a correlation analysis, long-term 
observations are essential. While LISA will not allow an opportunity for such 
a correlation analysis, some concept studies for the Big Bang Observer consider 
two systems such that the improvement related to the correlation analysis can 
be exploited. 

In Fig. 5, we concentrate on the frequency regime relevant for space-based 
detectors. In addition to detector noise and the expected level of amplitude 
related to inflationary gravitational waves, based on the energy scale of inflation, 
we also show the confusion related to foreground stochastic backgrounds related 
to merging compact binaries. These signals include cosmological neutron star- 
neutron star binaries (e.g., Schneider et ah 2001), white dwarf-white dwarf 
binaries (e.g.. Farmer &: Phinney 2003), neutron star-white dwarf binaries (e.g., 
Cooray 2004d), and the background related to single rotational neutron stars 
(Ferrari et ah 1999). In addition to these compact binaries, a cosmological 
background of hypothetical Pop III supernovae could potentially contribute a 
substantial background at frequencies corresponding to the Big Bang Observer 
(e.g., Buonanno et ah 2004). 

To jointly constrain the inflaton potential using observations at horizon- 
scale wavelengths with CMB and million-km scale with space-based interfer¬ 
ometers, one can extend the CMB-based technique discussed in Turner (1993). 
The detectable amplitude, and the spectral index, with space-based detectors 
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Figure 4. Here we show the predicted gravitational wave background 
due to inflation along with various current and planned detectors and 
known upper limits on a cosmological stochastic gravitational wave 
background. We follow Maggiore (2000) and take the correlated sen¬ 
sitivity to be approximately lO”'^ better than the uncorrelated sensi¬ 
tivity. The sensitivity of LISA is obtained from Shane Larson’s web¬ 
site, http://www.srl.caltech.edu/~shane/sensitivity/. The sensitivity 
for BBO correlated is obtained from Buonanno et al. (2004). We esti¬ 
mate the sensitivity for BBO by increasing the BBO curve by 4 orders 
of magnitude. The sensitivity for the ultimate DECIGO is estimated 
in Seto et al. (2001). The BBN constraint comes from the requirement 
that the total energy density in gravitational waves be less than the 
energy density in other forms of radiation at BBN, so that the abun¬ 
dances of the light elements remain in agreement with observations 
(Allen 1996). Only length scales less than the horizon size at BBO 
contribute and are constrained by this argument. Millisecond pulsar 
timing observations have constrained the energy density in a stochastic 
background to be below ~ 10“® around / ~ 10“® Hz. Curves corre¬ 
sponding to an inflationary produced background are shown in ma¬ 
genta for three different inflationary energy scales, (Turner, White 
&: Lidsey 1993). As stated in the text, COBE/WMAP constrain the 
energy scale of inflation to be below approximately, 2 x 10^® GeV. As 
is described in the text, due to confusion with weak leasing, a detec¬ 
tion of tensor perturbations in the polarization of the CMB can occur 
only if the energy scale of inflation is greater than ~ 2 X IQi® GeV. 
As of now, there is no planned gravitational wave observatory capable 
of detecting a stochastic background for energies at this lower energy 
bound. 


can be used to establish the inflaton potential amplitude and its first derivative. 
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Figure 5. Same as Fig. 4, but we show the background at frequencies 
related to the proposed Big Bang Observer mission. For comparison, 
we also show the expected stochastic backgrounds from foreground 
sources such as neutron star-neutron star binaries (e.g., Schneider et 
ah 2001), white dwarf-white dwarf binaries (e.g.. Farmer & Phinney 
2003), neutron star-white dwarf binaries (e.g., Cooray 2004d), and the 
background related to single rotational neutron stars (Ferrari et ah 
1999). For a reliable detection of the inflationary gravitational wave 
signal, the background related to neutron stars must be removed. This 
can be achieved given the low merger rate of these binaries and the 
rapid orbital evolution of the binary such that the gravitational wave 
frequency across this band changes substantially over the observational 
period. 


Combining this information with CMB, at roughly 30 orders of magnitude dif¬ 
ference in scales, should provide additional information to reconstruct the shape 
of the potential. While a priori assumed potential can easily be constrained 
based on the shape of the potential al, for a general reconstruction of the po- 
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tential shape, the simple approximations considered in the literature, based on 
power-law expansions of the potential, however, may not be appropriate given 
the large difference in scales. A detailed analysis on the extent to which in¬ 
flation models can be studies with joint constraints between CMB and direct 
gravitational-wave detections will be discussed elsewhere (Smith et ah 2005). 
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